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Abstract 

It has been known for some time that the Boltzmann weights of the chiral 
Potts model can be parametrized in terms of hyperelliptic functions, but as 
yet no such parametrization has been applied to the partition and correlation 
functions. Here we show that for N = 3 the function S(t p ) that occurs in the 
recent calculation of the order parameters can be expressed quite simply in 
terms of such a parametrization. 

1 Introduction 

There are a few two-dimensional models (and even fewer three-dimensional 
models) in equilibrium statistical mechanics that have been solved exactly. 
These are lattice models where spins crj are assigned to the sites i of a lattice 
(usually the square lattice). Each spin takes one of N possible values and spins 
<7i,Oj on adjacent sites i,j interact with a specified positive real Boltzmann 
weight function W((jj, crj). One wants to calculate the partition function (also 
called the sum-over-states) 

Z = E II W(a jlt r 3 ) , (1) 

<ij> 

where the sum is over all states of all the spins, and the product is over all 
edges of the lattice. 

If the number of sites is M, we expect the limit 

K= lim Z l ' M (2) 

to exist and to independent of the shape of the lattice, provided it becomes 
large in all directions: this is the "thermodynamic limit", and k is the expo- 
nential of the free energy per site. If 1, . . . ,m are sites fixed on the lattice 
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and the limit is taken so they become infinitely deep within it, then we also 
expect the average 

{f <r m )) = ^ _1 E II f(<n,---,*m)W(a j ,a j ) (3) 

<ij> 

to tend to a limit, for any given function / of these m spins. 

Because spins only interact with their neighbours, one can build up the 
lattice one row at a time, and associate a row-to-row "transfer matrix" with 
such an operation. 

To solve such a model, typically one shows that the Boltzmann weights W 
satisfy the star-triangle or "Yang-Baxter" relations [I]. These ensure certain 
commutation relations between the transfer matrices, and this is usually a 
first step towards calculating k. 

The next step is to calculate the order parameters, which are averages of 
certain functions of a single spin a\ deep within the lattice. This is a harder 
problem than calculating k. For instance, Onsager |2f)j calculated k for the 
square-lattice Ising model in 1944 , but it was not till 1949 that he announced 
at a conference his result for the order parameter (namely the spontaneous 
magnetization), and not till 1952 before a proof of the result was published 
by Yang [21] . 

However, since then the "corner transfer matrix" method has been devel- 
oped by Baxter [3], and the "broken rapidity line method" by Jimbo, Miwa 
and Nakayashiki ^J]- For many of the solved models (those with the "ra- 
pidity difference" property), these methods make the calculation of the order 
parameters comparitively straighforward. 

Even so, one model has proved challenging, namely the chiral Potts model. 
This is an iV-state model where W(aj,aj) depends only on the spin difference 
Gj — aj, mod N. The Boltzmann weights also depend on two parameters p, q, 
(known as "rapidities"), and on given positive real constants k, k' , related by 

k 2 + k' 2 = l . (4) 

The parameter k! plays the role of a temperature, being small at low temper- 
atures. For < k' < 1 the system displays spontaneous ferromagnetic order, 
becoming critical as k' — > 1. 

Its order parameters can be taken to be 

Mr = (lo™ 1 ) , (5) 

where u> = exp(2-7ri/iV) and r = 1, . . . , N — 1. It was shown in 1988 that its 
Boltzmann weights satisfy the star-triangle relation 012]) an d the partition 
function per site n was soon calculated |H] . 

The order parameters were another story. The model had developed from 
a one-dimensional quantum spin chain, which has the same order parameters. 
From series expansions it was conjectured pQ in 1989 that 

M r = F^"^ 2 . (6) 

Much effort was expended in the ensuing years (certainly by the author) in 
attempting to derive this result. It was not until 2005 that this was done 

PUCE]. 
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The method used was based on that of Jimbo et al [BJ] . In Figure Q we 
show the square lattice £, drawn diagonally, denoting the sites by circles and 
the edges by solid lines. We also show as dotted (or broken) lines the medial 
graph of C. Every edge of C is intersected by two dotted lines. With each 
dotted line we associate a rapidity variable (p, q, h or v). In general these 
variables may differ from dotted line to dotted line. They must be the same 
all along the line, except for the horizontal broken line immediately below the 
central spin a\ . We break this below c\ and assign a rapidity p to the left of 
the break, a rapidity q to the right. With these choices of rapidities, define 

F pq (r) = (uO . (7) 

In the thermodynamic limit, the star-triangle relations will ensure that 
F pq (r) is independent of the " background" rapidities v, h, because it allows 
us to move any of these dotted lines off to infinity. [H] However, the effect of 
the break is that we cannot move the broken line p, q away from <n, so F pq (r) 
will indeed depend on p and q. 

An important special case is when q = p. Then the p, q rapidity line is not 
in fact broken, so it can be removed to infinity and F pp (r) must be independent 
of p and equal to the order parameter A4 r defined by ® : 

Mr = F pp (r) . (8) 
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We also define 

G Pq (r) = F pq (r)/F pq (r-1) . (9) 

The author wrote down functional relations satisfied by G pq (r) in 
1998. They do not completely specify G pg (r), but must be supplemented 
by information on the analyticity properties of G pq {r). (Just as the relation 
f(z + 1) = f{z) only tells us that f(z) is periodic of period 1: however, if we 
can also show that f{z) is analytic and bounded in the domain < < 1, 
then it follows from Liouville's theorem that f(z) is a constant.) 

For N = 2 the chiral Potts model reduces to the Ising model and it is 
quite easy to find the needed analyticity information, to solve the functional 
relations and obtain the Onsager-Yang result M± = fc 1//4 . 

For N > 2 the problem is much harder. It was not until late 2004 that the 
author realised that it is not actually necessary to solve for the general function 
G pq (r). It is sufficient to do so for a special "superintegrable" case where q 
is related to p. The function then has quite simple analyticity properties 
and it quite easy to solve the relations (in fact one does not even need all 
the relations), to obtain G pq {r) for this case and to verify the 16-year old 
conjecture ©• For r = 1, . . . , N, the functions G pq {r) can all be expressed in 
terms of a single function S(t p ) which is defined below. 

Even so, it would still be interesting to understand G pq (r) more gener- 
ally. A fundamental difficulty is that for N > 2 the rapidities p and q are 
points on an algebraic curve of genus greater than 2, and there is no explicit 
parametrization of this curve in terms of single-valued functions of a single 
variable. (There is for N = 2: one can then parametrize in terms of Jacobi 
elliptic functions.) One can parametrize in terms of hyperelliptic functions [§], 
but these have N — 1 arguments that are related to one another. As yet they 
have not proved particularly useful, but one lives in hope. The function S(t p ) 
is a simple example of a thermodynamic property of the chiral Potts model, 
and has the simplifying feature that it depends on only one rapidity, rather 
than two. It is an interesting question whether it can be simply expressed in 
terms of these hyperelliptic functions. 

For N = 3 these hyperelliptic functions can be expressed in terms of 
ordinary Jacobi elliptic functions. One still has two related arguments (here 
termed z p and w p ), but some of the properties can be expressed as products 
of Jacobi functions, each with an argument z p or w p , or some combination 
thereof. A number of such results have been obtained. |TQ|, PP- 568, 569] 

There are two distinct ways of performing the hyperelliptic parametriza- 
tion. In |5J we used what we shall herein call the "original" parametriza- 
tion. What we report here is that for N = 3 the function S(t p ) can be 
expressed quite simply as a product of Jacobi functions of z p and w p , provided 
we use the second "alternative" parametrization. 

2 The function S(t p ) 

We can take a rapidity p to be a set of variables p = {x p , y p , /i p , t p } related to 
one another by 

tp = x p y p , x p + y p = k{\ + Xp i/p ) , 
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(10) 

kx p = 1 k I , ^Up == 1 ^ l^p 

There are various automorphisms or maps that take one set {x p , y p , n p ,t p } 
to another set satisfying the same relations (|1U() . Four that we shall use are: 

R-{x Rp ,y Rp ,fi Rp ,t Rp } = {y p ,ux p ,l/fip,ujtp} , 

S : {xsp,ysp,VSp,ts P } = {y p l ,x p l ,u- l/2 y p /(x P Li p ),t~ 1 } , 

V : {x V p,y V p,n V p,t V p} = {xp,ujy p , fi p ,ujtp} , (11) 

M : {x M p,yMp,HMp,t M p} = {x p ,y p ,ujfip,ujtp} . 

They satisfy 

RV~ 1 R = V , MRM = R , MSM = S , 

S 2 = V N = M N = 1 . (12) 
Let q be another rapidity set, related to p by q = Vp, i.e. 

x q = x p , y q = ojy p , fj, q = ^ . (13) 

We take fi p to be outside the unit circle, so 

\» p \ > 1 . (14) 

Then we can specify x p uniquely by requiring that 

- ?r/(2JV) < arg(xp) < ir/(2N) . (15) 

We regard x p ,y p ,/j,p as functions of t p . Then t p lies in a complex plane 
containing N branch cuts £>o, B±, . . . , Bn-i on the lines arg(i p ) = 0, 2tt/N, . . . , 
2tt(N — 1)/N, as indicated in Fig. [2J while x p lies in a near-circular region 
round the point x p = 1, as indicated schematically by the region IZq inside the 
dotted curve of Fig. |2J The variable y p can lie anywhere in the complex plane 
except in TZq and in N — 1 corresponding near-circular regions TZ±, . . . , 1Zn-i 
round the other branch cuts. With these choices, we say that p lies in the 
"domain" V. 

With these choices, we show in ^7] that 

G pq {r) = k^ + ^/ N2 S(t p ) , (16) 

for r = 1, . . . , N — 1, while 

G Pq (0) = Gp q (N) = k^' N2 Sitp) 1 ^ . (17) 

Hence G pq {\) ■ ■ ■ G pq (N) = 1, in agreement with the definition @. The func- 
tion S p = S(t p ) is given by 

logS(t p ) = ~ W2 ^k + — j o T - ¥ ^ log[A(60-tp]d* , (18) 
where 

A(0) = [(l-2k'cose + k' 2 )/k 2 ] 1/N . (19) 
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Figure 2: The cut t p -plane for N = 3, showing the three branch cuts B ,Bi, B 2 and 
the approximately circular region TZq in which x p lies when p G T>. 

From |T|2], particular properties are 

5(0) = 1 , 5(oo) = k- 2 ' N2 , 

S{t p )S{ujt p ) ■ ■ ■ S{uj n -\) = k~ l / N x p . (20) 

The function S{t p ) is single-valued, non-zero and analytic in the cut t p plane 
of Figure |2J but only the cut on the positive real axis is necessary: the other 
cuts can be removed for this function. If S ac (t p ) is the analytic continuation 
of S(t p ) across the branch cut B r , then 



S ac (t p ) = S(t p ) for r / 

= (y p /x p )S(t p ) for r = . 



(21) 



If we interchange p, q in eqn. 49 of ^7j, then apply the restriction (|13|) 
and use the relation RS = MVRSV together with eqn. 53 of J7]j we obtain 



GpJr)Gp,AN -r) = l , 



(22) 



where p' = V q' = RSVp. It follows that S(t p ) also has the symmetry 



S p Srsvp = S(tp)S(l/t p ) = k 2 I N . 



(23) 
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3 The Riemann sheets ( "domains" ) formed 
by analytic continuation 

We shall want to consider the analytic continuation of certain functions of t„ 
onto other Riemann sheets, i.e. beyond the domain D. We restrict attention to 
functions that are meromorphic and single- valued in the cut plane of Figure [3 
and similarly for their analytic continuations. Obvious examples are x p , y p and 
S(t p ). They are therefore meromorphic and single- valued on their Riemann 
surfaces, but we need to know what these surfaces are. 

We start by considering the most general such surface. As a first step, allow 
Hp to move from outside the unit circle to inside. Then t p will cross one of 
the N branch cuts Bi in Figure [21 moving onto another Riemann sheet, going 
back to its original value but now with y p in IZi . Since y p is thereby confined 
to the region near and surrounding u> 1 , we say that y p ~ u>' 1 . Conversely, by 
y p ~ uj % we mean that y p € TZi 

We say that p has moved into the domain T>i adjacent to T>. There are N 
such domains T>q,T>i, . . . , T>n-i- 

Now allow fi p to become larger than one, so t p again crosses one of the 
TV branch cuts. Again we require that t p returns to its original value. If it 
crosses Bi, then it moves back to the original domain T>. However, if it crosses 
another cut Bj then x p moves into TZj-i, and we say that p is now in domain 

Proceeding in this way, we build up a Cayley tree of domains. For instance, 
the domain T>ijk is a third neighbour of V, linked via the first neighbour T>i 
and the second-neighbour T>ij, as indicated in Figure El Here x p ~ 1 in T>, 
y p ~ w' in T>i, x p ~ uji in T>ij and y p ~ ui k in T>ijk- We reject moves that 
take p back to the domain immediately before the last, so j ^ and k 7^ i. 
We refer to the sequence {i,j, k, . . .} that define any domain as a route. We 
can think of it as a sequence of points, all with the same value of t p , on the 
successive Reimann sheets or domains. 

The domains T>, D{j, £>j^,... with an even number of indices, have x ~ ur, 
where £ is the last index. We refer to them as being of even parity and of type 
I. The domains T>i, T>ijk,... have y ~ ur and are of odd parity and type t. 



V T>i T>ij Vijk 



Figure 3: A sequence of adjacent domains V, Di, Vij,Vijk- 

The automorphism that takes a point p in T> to a point in X>j, respectively, 
is the mapping 

Ai = V 1 " 1 RV 1 . (24) 



If q = Ai p, then 



x q = uj l y p , y q = uj 1 x p , t q = t p . (25) 

Because of p2|). Aj + 7v = A,, so there are A" such automorphisms. 

We can use these maps to generate all the sheets in the full Cayley tree. 
Suppose we have a domain with route {i, j, k, . . .} and we apply the automor- 
phism A a to all points on the route. From (|25|) this will generate a new route 
{a, i — a,j + a, k — a, . . .}. For instance, if we apply the map A a to the route 
{m} from V to T> m , we obtain the route {a,m — a} to the domain D oljm - ct . 
Thus the map that takes T> to T>ij is AiAi + j. 

Iterating, we find that the map that takes T> to T>ijk ^ mn is 

AiAi+jAj+k ■ ■ ■ A m + n . (26) 

We must have 

4 = 1, (27) 

since applying the same map twice merely returns p to the previous domain. 

Let us refer to the general Riemann surface we have just described as Q. 
It consists of infinitely many Riemann sheets, each sheet corresponding to a 
site on a Cayley tree, adjacent sheets corresponding to adjacent points on the 
tree. A Cayley tree is a huge graph: it contains no circuits and is infinitely 
dimensional, needing infinitely many integers to specify all its sites. 

Any given function will have a Riemann surface that can be obtained from 
Q by identifying certain sites with one another, thereby creating circuits and 
usually reducing the graph to one of finite dimensionality. 

From ((25|) . the maps Aq, Ai, . . . , Ajy^i leave t p unchanged. We shall often 
find it helpful to regard t p as a fixed complex number, the same in all domains, 
and to consider the corresponding values of x p , y p (and the hyperelliptic vari- 
ables z p , w p ) in the various domains. To within factors of u, the variables x p 
and y p will be the same as those for T> in even domains, while they will be 
interchanged on odd domains. 

Analytic continuation of S(t p ) 

Now return to considering the function S(t p ). It is sometimes helpful to write 
this more explicitly as S{x p ,y p ). Then from (|21j) the map that takes S{t p ) 
from domain T> to T>i is 

q = Aip : S(x g ,y q ) = {y q /x q )~ &l S(x p , y p ) , (28) 

where x q ,y q are given by (|25|) and Si = 1 if i = 0, mod A^; otherwise 5i = 0. 
Note that x q ,y q are obtained by interchanging x p ,y p and multiplying them 
by powers of uj. 

For given t p , let So(t p ) be the value of S(t p ) in the central domain T>, given 
by the formula (|18|) . Iterating the mappings (|28|) from domain to domain, in 
any domain we must have 

S(x p ,y p ) = u a (y p /x p ) r S (t p ) , (29) 



S 



where a, r are integers. Note that in this equation x p , y p are the values for 
the domain being considered: they are not the corresponding initial values of 
the central domain T>. 

In particular, in the domain T)^ we obtain 

r = -Si + 5 i+j - 5 j+k . (30) 

4 The original hyperelliptic parameteriza- 
tion for N = 3 

Hereinafter we restrict our attention to the case N = 3 and use the hyperel- 
liptic parametrization and notation of previous papers. E3 EEU EE] We use 
only formulae that involve ordinary Jacobi elliptic (or similar) functions of 
one variable. 

Given k, k' , we define a "nome" x by 

oo /, _ 3n\ 12 

We regard x as a given constant, not the same as the rapidity variable x p . It 
is small at low temperatures (k' small), and increases to unity at criticality 
(k' = 1). We introduce two elliptic- type functions 

2w , ™ (l-ujx n - 1 z)(l-u 2 x n /z) 
- *V3 TT (1 - X 3 "" 2 /.) (1 - g^£) 

" n 1 = 1 1 (l-x3«-2 z)(1 _ x 3n-l /z) • W 

We then define two further variables z p , w p by 

t p = x p y p = ujh(z p ) = h(-l/w p ) = uj 2 h(-w p /z p ) , (34) 
These are the relations (27) of ^U|- The relations (32) of ^U] are also satisfied: 
x p 3 Vp /V 6 = H xz p/ w I) 3 = (^{-xZpWpf = cp(-xw p /z p ) 3 , (35) 
as are the relations (4.5), (4.6) of [TT] . in particular 



oo 

W 



r (1 - x 2n - l z/w){l - x 2n - x wjz){\ - x 6n - 5 zw)(l - x^^z^w' 1 ) 
W (1 - x 2n ~ 2 z/w)(l - x 2n w/z){\ - x 6n - 2 zw)(l - x 6 ™- 4 ^- 1 ^- 1 ) 

(36) 

writing z p , w p here simply as z,w. 

The z p , w p variables satisfy the automorphisms 

z Rp = xzp , z Sp = l/(xz p ) , Z Vp = -l/Wp , 

WRp = z p /w p , w Sp = l/(xw p ) , w V p = Zp/w p . (37) 
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The operation p — > Mp multiplies (z p w p ) 1 ^ by lo, but does not change z p , w p 
themselves. 

The variables z p , w p are of order unity when k', x are small, \jl p is of order 
1/k', and x p ~ 1. This is the low-temperature limiting case of p € D. It is 
convenient to define 



{z p , -1/Wp, -Wp/Zp} 



(38) 



The three automorphisms that leave t p unchanged, while taking T> to 
Dn,£>i,£>2, respectively, are 



-4,, 



V 2 R 



A, = RV 2 , A 2 = VRV 



If q = A{ p, then 



(39) 



(40) 



and u q = AiU p , where Aq,Ai,Ai are the three-by-three matrices 



An 



( 



\ 






X 1 





X 














1 



/ 









X 





1 





-1 









\ 



( I 



A, 



1 














X~ l 





X 






They satisfy the identities 



>Ai *Aj *A.i 



•^^■j '^^■'i 



(41) 



for all 

They permute the three elements z p ,—l/w p ,—w p /z p of u p and multiply 
them by powers of x, the product of the elements remaining unity. Let z p , w p 
be the values of z p , w p on the central sheet T>. Then it follows that on any 
sheet, for the same common value of t p , 



x a 



p ' 



i". 



x n f3 p 



(42) 



where {a p , —l/(3 p , —f5 p /a p } is a permutation of {z p , —l/w p , — w p /z p }. 

Repeated applications of the three automorphisms will therefore generate 
a two-dimensional set of permutations and multiplications of the elements of 
u p . Each member of the set corresponds to a site on the honeycomb lattice of 
Figure HI Adjacent Riemann sheets correspond to adjacent sites of the lattice. 
Sheets of even parity correspond to sites represented by circles, those of odd 
parity are represented by squares, if i is the integer inside the circle or square, 
then for even sites x p ~ while on odd sites y p ~ oj % . The numbers shown 
in brackets alongside each site are the integers m,n of (I42|) . 

Thus for the functions z p and w p of t p , the graph Q of the Riemann surface 
reduces to this two-dimensional honeycomb lattice. 

Note that the sites X, Y, Z in the figure are third neighbours of the central 
site V, and each can be reached from T> in two three-step ways. For instance, 
Y is both £>n2i and X^n- 1 



1 Note that for T>ijk we here take the intermediate site j to be represented Figure 0] by 
integer — j, mod 3. This is changed in the next section to +j. 
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Thus Y is obtainable from V by the maps A0A2A0 and A2A0A2. From 
(|4ip these are the same, so we can identify the two sheets as one, represented 
by Y. 

Similarly, X corresponds to A1A2A1 = A2A1A2 and Z to A0A1A0 = 
j4i.Ao.Ai. This is why Q reduces to the honeycomb lattice. 

However, the automorphisms (|28|) of the function S(t p ) do not in general 
satisfy (|41|). On its Riemann sheets P02I) ^211 we find from (|3fl|) that r = —2 
and 1. Thus from the analytic continuation of S(t p ) is (y p /x p )~ 2 So(t p ) 
and (y p /x p )So(t p ) on each sheet, respectively (ignoring factors of Thus 
the result for Y depends on the route taken to it. (The same is true of Z, but 
not for X.) 

Hence the values of z p (and w p ) are the same on D021 and D21I1 but 
S(t p ) is different. It follows that z p , w p do not uniquely determine S(t p ). 
Hence neither S(t p ) nor S(t p ) 3 is a single- valued function of these hyperelliptic 
variables z p ,w p : one must look elsewhere for such a parametrization. 



X (2,2) 



U (-2,0) 




2 (-3,-2) 



1 (-2,-3) 



Figure 4: The honeycomb lattice formed by the hyperelliptic variables z, w in the alternative 
parametrization. Circles (squares) denote sites of even (odd) parity. 
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5 The alternative hyperelliptic parametriza- 
tion for N = 3 



There is another way of parametrizing k,x p ,y p , n p ,t p so that the nome x is 
small when k' is small. It can be obtained from the original parametrization 
of the previous section by a simple mapping, and we do this in the Appendix. 
We take the results (|A2|) . (|A3|> . (|A4|I . (|A3|) therein and drop the hats on 
x p , y p , . . . , V, M to obtain 

-*' 2 =^n(^j • <43) 

coh(zp) = h(—l/wp) = uj 2 h(—w p / z p ) , (44) 
0(xz p /wp) 3 = 0(-xZpWJp) 3 = ^)(-xu;p/Zp) 3 . (45) 



ZRp = -xw p , z Sp = -l/(xw p ) , Z V p = -l/Wp , 

WRp = W p /Zp , W S p = -1/(XZ P ) , W V p = Z p /Wp . (46) 

Again zm p = z p , WMp = w p . These equations replace (J3*T|) . (j35j) . (|3T|) of 
section 4. The functions /j(z), 0(^) remain defined by (|H2*|) and (|3*H|) . and the 
relation ()36|) remains satisfied. 

We now regard x p ,y p , . . . ,V, M as being the same variables and automor- 
phisms as those above, satisfying (fTU)) . (fTTj) and (|12|) . Then the hyperelliptic 
variables e different from those of section 4. If one takes p £ T> 

and expands the functions in powers of the low-temperature variable k', then 
to leading order the same as the z p , w p of section 4, being of order 

unity and satisfying w p = z p + 1; x is negated and is of order /c'. 

The three automorphisms that leave t p unchanged while taking P to Do, 
T>i, T>2 are again given by (f3H)) and (|IUJ). Using the rules we find that 
Zp,w p transform according the rules 

q = A lP : z g = x 2 ^- 3<5 «/^ , ^ = x'-'/wp , (47) 

for i = 0, 1, 2, writing <5j as 5(i). If z^, are the values of z p , w p on the central 
sheet T>, then it follows that on any Riemann sheet the analytic continuations 
of Zp, Wp (for a given value of tp) are 

x m (z a p ) ±l , x n {w§ ±l , (48) 

choosing the upper (lower) signs on sheets of even (odd) parity. Here to, n are 
integers satisfying 

to + n = (mod 3) on even sheets , 

= 1 (mod 3) on odd sheets . (49) 

The Riemann surface for z p , w p therefore corresponds to a two-dimensional 
graph Q, each site of Q being specified by the two integers m, n. 
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In fact this Q is the honeycomb lattice shown in Figure 0J but we must 
interpret it slightly differently from what we did in section 4. Adjacent sites 
still correspond to adjacent Riemann sheets, and the parities of the sites are 
shown as in section 4 by circles and squares. Again, on odd sites y p ~ w ! , 
where i is the number shown in the Figure. However, on even sites we now 
take x p ~ lo 1 . For each site, the bracketed integers shown in Figure |1] are now 
the integers (m,n) of (fl8|) . 

As in the previous section, the sites X, Y, Z in the figure are third neigh- 
bours of the central site P, and each can be reached from T> in two three-step 
ways. There is the difference that the j of T>^ k is now the number inside 
the corresponding circle of Figure 01 so in this section X is both D210 an d 
^120) s° from (|2fij) and (|27j) it is obtainable from T> by the maps A2AQA1 and 
A\AqA<i- The reason we can identify D210 with X'120 is that ()47|) implies that 
AiA§A\ = A\A§A2- More generally, it implies that 

A t AjA k = A k AjAi (50) 

for all permutations i, j, k of 0, 1, 2. This means that X, Y, Z each corresponds 
to a single Riemann sheet rather than two, and is the reason that Q reduces 
(for the functions z p ,w p of t p ) from the full Cayley tree to the honeycomb 
lattice. 



Properties of z p , w p , S(t p ) 

Within the central domain T> there is some circle of non-zero radius, centre 
the origin, such that none of z p ,w p , l/z p ,l/w p , Zp/w p ,Wp/z p can lie within 
the circle (for x small the radius is of order x 1 / 2 ). Two special values of p 
that lie within T> are 

p(l) : z p = -l/w p = -Wpjzp = J 1 , y p = 

x p = k 1 ^ , S(t p ) = 1 , (51) 

and 

p(2) : z p = -l/wp = -Wp/zp = lo , y p = oo 

using eqns. ((TU)) and (|2"U)) . 

Remembering that x p ~ 1, one can verify that the conditions (j.36j) . (|44|) 
are satisfied. These are the only two points within T> where y p /x p has a zero 
or pole. 

Any point where one of z p , —l/w p , —w p /z p is equal to u> 2 x m or ujx™ 1 , for 
non-zero integer m, necessarily lies outside T> (i.e. on another Riemann sheet). 
It follows that the functions G(z p ), . . . , h{—Wp/z p ) defined and used in the next 
section have no zeros or poles for p GP: they are finite and non-zero therein. 
The same is true of x p , S(t p ) and the function S(z p , w p ) defined below. 

Now consider the maps Ai for the function S(t p ), as given in (|28|) - (|3U|) . 
The mapping A^AjA^ simply multiplies S(t p ) by a power of u; and by (y p /x p ) r , 
where r = —5i + 5j — 5 k - Obviously this exponent r is unchanged by inter- 
changing i with k, so is the same for both sides of (l50l) . 
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The powers of u> are not necessarily the same, but we can avoid this diffi- 
culty by simply working with S(t p ) 3 instead of S(t p ). The automorphisms 
for S(t p ) 3 do satisfy (jHOjl and S(t p ) 3 is the same on the two sheets X>2io> £>i20) 
so again X reduces to a single sheet. Similarly, so do Y and Z. 

As a result, S(t p ) 3 is uniquely determined if we know z p and w p . In fact 
on any Riemann sheet (m, n) the analytic continuation of S{t p ) 3 is 

(y P /x p ) 3r S (t p f , (53) 

where 

r = (m + n)/3 on even sheets , 

= (m + n— l)/3 on odd sheets . (54) 

The x p , y p in (|53j) are those of the sheet under consideration. 

It therefore makes sense to look for a single-valued meromorphic function 
of z p and w p that is equal to S(t p ) 3 . We can write S(t p ) itself as S(z p ,w p ) 
provided we accept that it is three-valued in the rather trivial way that its 
values differ by factors of lj. Then from the automorphisms ()28|) and (|47|) . 
this function must satisfy the three relations 

S(z p ,Wp) = (xp/y^Stx" 1 z^ 1 , x" 1 ^" 1 ) =5(iz~ 1 ,w p " 1 ) = S{z~ x , x w" 1 ) 

(55) 

for all Zp, Wp. Also, using JTTJ) and (JS)), the relation ((201) becomes 

S(zp,Wp)S(-l/wp,Zp/wp)S(-Wp/zp,-l/z p ) = k~ 1/3 x p . (56) 

The function x~ l has no zeros or poles on even sheets, whereas on odd 
sheets (where y p ~ oj % ) it has the same zeros and poles as y p /x p . it is therefore 
useful to work not with S(z p ,w p ) 3 , but the function 

S(zp,w p ) = x~ l S{zp,Wp) 3 , (57) 

since from (|5T?|) and (|5*I|) this has the same poles and zeros as (y p /xp) m+n on 
all Riemann sheets (m,n), even and odd. From we can write (y p /x p ) as 
a function of either z p ,w p or z p /w p : this suggests that it may be possible to 
write S(z p , w p ) as a product of functions of these individual variables. We do 
this in the next section. 



6 S(tp) as a function of z p , w t 



Define the functions 



*"(*) = n i 



(i - xHy 



1J ; (1 - xiz~ l )i 

G(s) = F(luz)/F(lu 2 z) , (58) 
so F(z) = 1/F(z~ 1 ), G(z) = G{z~ l ) and 

G{z)/G{xz) = ujh(z) (59) 
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for all complex number z. 
Consider the product 

P = G(z p ) a G(-l/w p fG(-w p /z p y (60) 

for arbitrary integers a, (3, 7. As a first step, we ask if we can choose a, /3, 7 
so that P has the same poles and zeros as S(z p , w p ), i.e. as (y P /x p ) m+n on all 
sheets (m,n). 

The integers m,n specify a sheet and are defined in (|48|). From (|44|1. on 
sheet (m, n) the function y p /x p has a simple zero when z p = x m uj 2 , w p = 
—x n tv, z p /w p = —LV m ~ n uj. Any one of these equalities implies the other two, 
so at this point S(z p ,w p ) has a zero of order m + n. On the other hand, the 
three factors of P have zeros of order —ma, n(3, (m — 71)7, respectively. Thus 
we require m + n = —ma + n(3 + (m — 71)7, for all allowed integers m, n and 
fixed values of a, /?, 7. This will be so iff 

a = 7-1 , /3 = 7 + l . (61) 

There are also possible zeros and poles when z p = x m uj, w p = —x n uj 2 , 
Zp/vjp = —uj m ~ n L0 2 , but the only difference from the above is that all the 
orders are negated, so again we obtain the conditions (|61jl. 

Neither S(z p ,w p ) not P has any other other zeros or poles, and S(z p ,w p ) 
is independent of the integer 7, which is still arbitrary. Substituting l|61|) into 
(|6U|) . this implies that the functions 

S(zp,Wp)G(zp)/G(-l/w p ) , G{zp)G{-l/wp)G{-Wp/z p ) 

have no zeros or poles anywhere on the Riemann surface. Using (|Blj) and (|57j) . 
they have values A;" 1 / 3 , F(uj)~ 3 at p = p(l), so if they were constants it would 
follow that 

S(zp,w p ) = k- l l*G{-l/wp)/G(zp) , (62) 

G{zp)G{-l/wp)G(-Wp/zp) = 1/F(lu) 3 . (63) 

We can prove that these relations are indeed true by using Liouville's 
theorem for a single Riemann sheet. Let T(z p , w p ) be the ratio of the RHS of 
(fHTjl to the RHS of (H2J). Then from and - it follows that 

T(z p ,w p ) =T(x~ 1 z p 1 ,x~ 1 w~ 1 ) = T(xz~ 1 ,w p 1 ) = T(z~ 1 ,xw~ 1 ) . (64) 

Thus the function T is unchanged by the three automorphisms Aq, A\, Ai- It 
is therefore a single- valued function of the variable t p , without the branch cuts 
of Figure [21 It has no zeros or poles in T>, so it has no zeros or poles in the 
complex t p plane, including the point at infinity. By Liouville's theorem it is 
therefore a constant. It is unity when z p = to 2 or lo, i.e. when t p = or 00, 
so it is one. This proves the identity ()62|) . The identity 1)63(1 can be proved 
similarly by taking T(z p , w p ) to be the ratio of of the LHS of l|63|) to the RHS. 
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Two more identities 



We originally tried a much more general ansatz for the product P, allowing 
for factors such as 1 — x° z raised to a power linear in j and mod (j, 3). As 
a result we discovered yet two more identities satisfied by z p ,w p . Define the 
function 

h( Z ) ■ n a-^-^xi-* 3 "' 1 ^) f65 x 

then we find that 

h(z p )h(-w p ) = h(-l/w p )h(-z p /wp) = h(-w p /z p ) h(l/z p ) . (66) 
We can prove these identities in a similar way. First note that 
h(x/z) = l/h(z) , h{z/x) =uj 2 /h(x- 1 z- 1 ) 

and 

h(z) = h{z/x)h{z)h{l/z) 

for all z. We can use this last formula and (|44|) to eliminate ratios such as 
h{— /h(x~ l z) in favour of h functions whose arguments do not contain 
x as a factor. Using this fact and applying the automorphisms Aq, A\,Ai to 
the ratios of the expressions in ()66|) . we find that the automorphisms merely 
permute these ratios. If we write the three expressions as J\, J2, J3 and form, 
for arbitrary a, 

T(z p ,w p ) = (a-J2/Ji)(a-J 3 /J2)(a-Ji/J 3 )(a-J 1 /J 2 )(a-J2/J3)(o!-J 3 /J 1 ) 

then this T(z p ,w p ) is unchanged by the automorphisms, so is a single-valued 
function of t p . It has no zeros or poles in T>, so by Liouville's theorem it is a 
constant. At p = p(l) or p(2), the Jj's are equal, so for all z p ,w p 

T(z p ,w p ) = (a - l) 6 . 

It follows that J\ = J2 = J3 for all t p , which establishes the identities (|66|) , 

For an arbitrarly chosen numerical value of z p , working to 32 digits of 
accuracy, we have successfully checked the identities (|62|) , (|63|) , (|66|) to twenty 
terms in an expansion in powers of x. 

From (|TT|) and (|l6*)) . the map p — > RSVp takes x p ,t p , z p ,w p to l/x p , l/t p , 
—w p , —z p . From (fHTj) and (|6^|) . noting that G(—l/w p ) = G(—w p ), it is ap- 
parent that 

S P S R sv P = S(z p ,w p )S(-w p ,-z p ) = k~ 2/9 , 
in agreement with ((2 



7 Summary 

For the case when N = 3 and q is related to p by (|13|) . the generalised 
order parameter function G pq (r) is given by (|Td]) . where x~ 1 5(i p ) 3 can be 
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simply expressed by (|57|) and (|62|) as a product of functions of the hyper elliptic 
variables z p ,w p of section 5. This is the first time that a thermodynamic 
property of the N > 2 chiral Potts model has been so expressed. (As distinct 
from algebraic functions of the Boltzmann weights, such as the function of f pq 

of H3|.) 

The functions F(z), G(z) in (|58(1 are infinite products similar to elliptic 
functions, except that factors such as 1 — x 3 z are raised to the power j. Such 
extensions of elliptic functions occur in the free energies of other solvable 
models, notably the Ising model. Good examples are the equations (B.10), 
(B.17) and (B.18) of 

We emphasize that (|57j) and (|62j) are in terms of the alternative hyperel- 
liptic parametrization of section 5 herein. Papers [0] to ^H] are in terms of 
the original hyperelliptic parametrization of section 4. 

It is still an interesting question whether G pq (r) can be simply expressed as 
a function of such variables for arbitrary p, q. The result of this paper implies 
that one must use the hyperelliptic parametrization of section 5, rather than 
that of section 4. There is a difficulty with this: if we write R 6 p, R e q as p', q', 
the relations (|46j) imply that z p i = z p and w p > = w p , whereas it is not true 
that G p ' iq (r) = G pq (r) or that G p ^ q i{r) = G pq (r). This means that G pq (r) 
cannot be a single-valued function of z p ,w p , z q ,w q . However, the function 
L pq( r ) = G pq (r)G Rq)Rp (r) of [T2] is unchanged by p -> p', and by q -> q', so 
may be so expressible. 



Appendix 

Here we show how the alternative hyperelliptic parametrization of section 5 
can be obtained from the original parametrization of section 4 by a simple 
mapping. 

Let k,k' ,x p ,y p , /j, p ,t p be the variables of section 4. Define new variables 
, k , x p ^ y P ] l^pi t p so that 

k = k , k = lk /k , x p = l/x p , y p = y p , 

fi p — e ^ fcpjj'p i yp/ %p — tp ■ (-^d) 

Leave x,z p ,w p and the functions h(z), 4>(z) unchanged. 

Then the relations ©, (JEJ) remain satisfied if we replace k, k' , x p , y p , t p 
therein by k,k' ,x p ,y p , p, p ,t p . The relations (|HT|) . (|H4jl . (jHSj) become 

y p /x p = ujh(z p ) = h(-l/w p ) = uj 2 h(-w p /z p ,x) , (A3) 
- x p 3 y p /l p 6 = 4>{xz p /w 2 p f = 4>(-xz p w p f = (j)(-xw p /z p f . (A4) 

Define automorphisms R, S, V, M by (jllj) with x p , y p , . . . , M p replaced by 
x p ,y p , . . .,M P . Then 

r = vs , 5 = M~ 1 i/~ 1 i? , y = y , M = M 
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and the relations (|12|) remain satisfied if R, S, V, M therein are replaced by 
R, S, V, M. From (j37|) it follows that 

Z Rp = ~ XW P ' Z Sp = - 1 I{XW P ) , Zy p = -1/Wp , 

W Rp = W p/ Z P ' W Sp = ~ 1 /( XZ p) > W Vp = Z p/ W P ■ ( A5 ) 

Now drop the hats on k, k , x p , . . . , M to obtain the equations (j%3*|) - 
of section 5. Equations (@J), (|10j) . remain true. 
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